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Abstract 

We prove a uniqueness result for limit cycles of the second order ODE 
x + i(f>(x, x) + g(x) — 0. Under mild additional conditions, we show that 
such a limit cycle attracts every non-constant solution. As a special case, 
we prove limit cycle's uniqueness for an ODE studied in [5] as a model 
of pedestrians' walk. This paper is an extension to equations with a non- 
linear g(x) of the results presented in [?]• 

Keywords: Uniqueness, limit cycle, second order ODE's, star-shaped 
function, Conti-Filippov transformation. 



1 Introduction 

The simplest non-linear continuous dynamical systems originate from the study 
of planar differential systems, 

x = P(x,y), y = Q(x,y), P,Q e ^(ia 2 ,^ 2 ). (1) 

Special cases of such systems are Lotka-Volterra ones, and systems equivalent 
to Lienard equations, 

x = y-F(x), y = -g(x), (2) 

or to Rayleigh equations 

i = y, y = -g(x) - f(y). (3) 
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All of them arise as mathematical models of biological, physical, engineering 
systems [3] . The study of the dynamics of ([I]) strongly depends on the existence 
and stability properties of special solutions such as equilibrium points and non- 
constant isolated periodic solutions. In particular, if an attracting non-constant 
periodic solution exists, then it dominates the dynamics of |l]) in an open, 
connected subset of the plane, its region of attraction. Studying the number 
and location of isolated periodic solutions, usually called limit cycles, is by no 
means a trivial question, as shown by the resistance of Hilbert XVI problem (see 
[9], problem 13). In some cases such a region of attraction can even extend to 
cover the whole plane, with the unique exception of an equilibrium point. In such 
a case the limit cycle is unique and dominates the system's dynamics, as in [J]. 
Uniqueness theorems for limit cycles have been extensively studied (see [2] , [T3] , 
|14j . for recent results and extensive bibliographies). Limit cycle's uniqueness 
is a relevant feature even in discrete time systems, which are often related to 
continuous time systems [10]. Sometimes, suitable symmetry conditions have 
been used, in order to simplify the study of such systems. In particular Z 2 
symmetry, that is orbital symmetry with respect to one axis, has proved to be 
useful in approaching similar problems |12j . 

Most of the results obtained for continuous time dynamical systems in the 
plane are concerned with the classical Lienard system ([2]) and its generalizations, 
such as 



x = £(x) 



<p{y) - F(x) 



y = -((y)g(x). (4) 



Such a class of systems also contain Lotka-Volterra systems and systems equiv- 
alent to Rayleigh equation Q as special cases. 

Even if the systems Q reach a high level of generality, compared to Van der 
Pol system, 

Even if the systems Q reach a high level of generality, compared to Van der 
Pol system, 



x = y-e y— - x\ , y=-x. 

the first one to be investigated in relation to existence and uniqueness of limit 
cycles, an evident limitation is given by the fact that the variables x and y 
appear separately, so that mixed terms are products of single- variable functions. 
Since models displaying a different combination of variables do exist, different 
methods are desirable, in particular in absence of symmetry conditions. 
A recent result [5] is concerned with systems equivalent to 

N 

x + ^2f 2k+1 (x)x 2k+1 +x = 0, (5) 

with / 2 fc+i(x) increasing for x > 0, decreasing for x < 0, k = 0, . . . , N. On the 
other hand, there exist classes of second order models which are not covered by 
previous results. This is the case of a model developed in [5] to describe the 
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pedestrian's walk, which leads to the equation 

x + ex(x 2 + xx + x 2 - 1) + x = 0, e > 0. (6) 

Such an equation can be considered as a special case of a more general class of 
equations, 

x + x(f>(x, x) + g{x) = 0. (7) 

In this paper we study the class |7]), assuming <p{x,x) to have strictly star- 
shaped level sets and xg{x) > for x =^ 0. We prove a uniqueness result for 
limit cycles, and, under suitable additional assumptions, we show that a limit 
cycle exists and attracts every non-constant solution. Since ex(x 2 + xx + x 2 — 1) 
has strictly star-shaped level sets, the model introduced in [5] has a unique limit 
cycle, attracting every non-constant solution. 

The result we present here is as well applicable to several equations of Lieanrd 
and Raylcigh type, in particular when they have a non-linear g(x). 

This paper is organized as follows. In section 1 we study the equation (m), 
assuming g{x) to be linear. We first prove the uniqueness theorem. The main 
tools applied here is a uniqueness result proved in . Then we introduce some 
mild additional hypotheses on the sign of <p(x,x), under which the unique limit 
cycle attracts every non-constant solution. Then, in section 2, we assume g{x) 
to be non-linear. We reduce the study of such a case to that of the linear g(x), 
by means of Conti-Filippov transformation |S]. The structure of section 2 is 
very similar to that of section 1, the main difference being the derivation of the 
condition on (f>(x,x) which implies the strict star-shapedness property for the 
transformed system. 



2 Linear g(x) 

Let fl C IR 2 be a star-shaped set. We denote partial derivatives by subscripts, i. 
e. (f> x is the derivative of 4> w. r. to x, etc.. We say that a function <f> G C 1 (il,IR) 
is star-shaped if (x, y) ■ V0 = x<p x + y4> y does not change sign. We say that <fr 
is strictly star-shaped if (a;, y) ■ V</> 7^ 0, except at the origin O = (0, 0). We say 
that j(t) is positively bounded if the semi-orbit 7+ = {"f(t), t > 0} is contained 
in a bounded set. Similarly for the negative boundedness. We say that an orbit 
is an open unbounded orbit if it is both positively and negatively unbounded. 
We say that a set X is invariant if every orbit starting at a point of X is entirely 
contained in X. For other definitions related to dynamical systems, we refere 
to pQ We call ray a half-line having origin at the point (0, 0). 

In this section we are concerned with the equation 

x + x<p(x, x) + kx = 0, k e IR, k > 0. (8) 

Without loss of generality, possibly performing a time rescaling, we may restrict 
to the case k = 1. Let us consider a system equivalent to the equation pi, for 
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k = 1, 

x = y y = -x - ycj>{x,y). (9) 

We denote by j(t , x* , y* ) the unique solution of the system ^ such that 7(0, x* , y* ) 
(x* 7 y*). We hrst consider a sufficient condition for limit cycles' uniqueness. We 
set 

A(x, y) = yx - xy = y 2 + x 2 + xycj>(x, y). 

The sign of A(x, y) is opposite to that of the angular speed of the solutions of 
([9| . Our uniqueness result comes from theorem 2 of [6] , in the form of corollary 
6. 

Theorem 1. Let cj> € (1R 2 ,IR 2 ) be a strictly star-shaped function. Then has 
at most one limit cycle. 

Proof. Without loss of generality, one may assume that, for (x,y) ^ (0,0), 
xcj) x + y<j) y > 0. The proof can be performed analogously for the opposite 
inequality. 

We claim that VA(x,y) does not vanish on the set Aq = {(x,y) : A(x,y) = 
0} \ {(0, 0)}. In fact, VA and A vanish simultaneously at (x,y) if and only if 

2x + ycj> + xy4> x = 

2y + xcf) + xycpy = . (10) 
x 2 + y 2 + xycf) = 

Multiplying the first equation by y, the second one by x and re-ordering terms 
yields 

xycj) = -x 2 (2 + ycf) x ) 
xycp = -y 2 (2 + xcj) y ) . 
xycj) — —x 2 — y 2 

Multiplying the third equation by 2 and summing with the first two equations 
yields 

xy(xcj) x + ycj) y ) = 0. 

Since, by hypothesis, xcj> x +ycj> y ^ except at O, one has xy = 0. If x = 0, then 



by the third equation in ( 10 ) one has y — 0. Similarly, if y = 0. 

This shows that, at every point, Aq is locally a graph. Additionally, every 
ray {{tcos9^tsm9),t > 0}, meets Aq at most at a point. In fact, for xy 7^ 0, 
one has, 

A(£cos0,tsin0) = <^=^ <j)U cos Q,tsvn6) = 1 . 

cos 6 sin 6 

The condition x<p x + ycp y > implies that cj> is an increasing function of t on 
every ray. Hence on every ray not contained in an axis there exists at most one 
t such that cj){t cos 6, tsin 9) = cos g 1 ^ g ■ As for xy = 0, A vanishes only at O. 
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Moreover, working as above, one can show that Aq has at a single point in 
common with the axes, O. 

The radial derivative A r of A is given by 

A r = xAx + r VAy =\ (2A + xy{x4> x + ■ (11) 

Let (x*,y*) be a point of the first orthant, i. e. x* > 0,y* > 0. If A(x*,y*) > Q, 
then A r > at and at every point (rx*,ry*) with r > 1, hence A is 

strictly increasing on the half-line {(rx*,ry*) : r > 1}. Now, let (x*,y*) be 
a point of the second orthant, i. e. x* > 0,y* < 0. If A(x*,y*) < 0, then 
A r < at (x*,y*) and at every point (rx*,ry*) with r > 1, hence A is strictly 
decreasing on the half-line {(rx*,ry*) : r > 1}. The same argument allows to 
prove that in the third orthant A behaves as in the first one, and in the fourth 
orthant A behaves as in the second one. 

Assume, by absurd, two distinct limit cycles to exist. The system Q has a 
unique critical point, hence they are concentric. Let [i\ be the inner one, /12 be 
the external one. Let D be the annular region bounded by fi± and We claim 
that A(x, y) > in D. We prove it by proving that, for every orbit 7 contained 
in D, A(j(t)) > 0. Let us observe that every orbit in D has to meet every semi- 
axis, otherwise its positive limit set would contain a critical point different from 
O. On every semi-axis one has A(x, y) > 0. Assume first, by absurd, A(7(i)) 
to change sign. Then there exist ti < t 2 such that 74(7(^1)) > 0, ^(7(^2)) < 0, 
and and r ){t{) 1 i — 1,2 are on the same ray. Assume "/(U), i = 1, 2 to be in the 
first orthant. Two cases can occur: either |7(ii)| < 17(^2)1 or |7(ii)| > |t(^2) | - 
The former, |7(^i)| < | T (*2 ) | , contradicts the fact that A is radially increasing 
in the first orthant, hence one has |7(ii)| > | t(*2 ) I - The orbit 7(^1) crosses the 
segment S = {rj(ti), < r < 1}, going towards the positive y-semi-axis. Let G 
be the region bounded by the positive y-semi-axis, the ray {r^(ti),r > 0} and 
the portions of fix, /12 meeting the y-axis and such a ray. The orbit 7 cannot 
remain in G, since in that case G would contain a critical point different from 
O. Also, 7 cannot leave G crossing the positive j/-semi-axis, because A(x, y) > 
on such an axis. Hence 7 leaves G passing again through the segment S. That 
implies the existence of i 3 > i 2 , such that 7(^3) lies on the ray {r^{t\), < r}. 
Again, one cannot have |t(*3) I < 17(^)1, since A(7(£ 3 )) > implies A increasing 
on the half-line rj(t 3 )),r > 1, hence one has 17(^3)! > |7(*2)|- Also, one cannot 
have |7(i3)| < |7(ti)|, otherwise 7 would enter a positively invariant region, 
bounded by the curve 7(i), for ti < i < t^, and by the segment with extrema 
7(^1), 7^3), hence there would exist a critical point different from O. As a 
consequence, one has |7(£s)| > |7(*i)|- Since A(x,y) > on the segment joining 
7(ii) and 7(^3), such a segment, with the portion of orbit joining 7(^1) and 7(^3) 
bounds a region which is negatively invariant for (j9|, hence contains a critical 
point different from O, contradiction. 

This argument may be adapted to treat also the case of a ray in the second 
orthant, replacing the positive y-semi-axis with the positive x-semi-axis, and 
reversing the relative positions of the points 7(ti), 7^2), 7(^3)- In the other 
orthants one repeats the arguments of the first and second orthants, respectively. 
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Finally, assume that A(j(t)) = at some point (x* , y*). If (x* , y*) is interior 
to D, then on the ray {(rx* ,ry*),r > 0}, there exist points interior to D with 
A(rx* ,ry*) < 0. Then we can apply the above argument to the orbits starting 
at such points. If (x*,y*) belongs to the boundary of D, then it is on /zi or 
on /i 2 . Assume (x* ,y*) — fj,i(t*), for some t* (the argument works similarly on 
^2). Since A(ni(t*)) = 0, fi% is tangent to the ray {{rx* , ry*), r > 0}. On the 

other hand, A r = — 7^ at (x*,y*), hence Aq is neither tangent 

r 

to the ray {(rx* ,ry*),r > 0}, nor to fi\ at fii(t*). This implies that Aq and fi\ 
are transversal at (x* , y*), so that a portion of A enters D. Since A separates 
points where A > from points where A < 0, also in this case there exist points 
interior to D with A < 0. 

Now we can restrict to the annular region D and divide the vector field of 
([9| by A(x, y), as in corollary 6 in [6]. In order to appy such a corollary, one 
has to compute the expression 

v = P (xQ x + yQ y ) - Q (xP x + yP y ) , 

where P and Q are the components of the considered vector field. For system 
([9]), one has 

v = y (-x - xy(f> x -y<j)- y 2 (t>y) - {-x -y<fr)y = 

-y 2 (x<t> x + y<t> y ) < 0. 

The function v vanishes only for y = 0. For both cycles one has: 



u(fa(t))dt < 0, i = 1,2, 

where Tj is the period of /Uj, i = 1,2. Hence both cycles, by theorem 1 in [5], 
are attractive. Let A\ be the region of attraction of fii. A\ is bounded, because 
it is enclosed by ^2, which is not attracted to fj,\. The external component of 
Ai's boundary is itself a cycle /Z3, because ([£]) has just one critical point at O. 
Again, 

p(fj, 3 (t))dt < 0, 







hence /13 is attractive, too. This contradicts the fact that the solutions of (|9| 
starting from its inner side are attracted to Hence the system ([9]) can have 
at most a single limit cycle. £ 

The angular velocity of the solutions need not be negative at every point 
of the plane. In fact, even a simple sysytem as that one studied in (|6| has 
negative angular velocity only in a proper subset of the plane. In Figure 1 we 
have plotted some orbits of the equation [6] tending at the limit cycle, together 
with the two components of the curve Aq. The orbits cross Aq at the points 
where their angular velocity changes sign. 
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In the example of Figure 1 we have chosen e = 1. In general, the system 
studied in jo] has just one limit cycle, for e > 0. In fact, in this case one has 

v = x(f) x + y<fi y = 2ey 2 (x 2 + xy + y 2 ) > for y ^ 0. 

It should be noted that even if the proof is essentially based on a stability 
argument, the divergence cannot be used in order to replace the function v. In 
fact, the divergence of system ^ is 

div (y, ~x - y<f)(x, y)j =-(/>- ycf> y , 

which does not have constant sign, under our assumptions. Moreover, the di- 
vergence cannot have constant sign in presence of a repelling critical point and 
an attracting cycle. 

Now we care about the existence of limit cycles. Let us denote by D r the 
disk {(x, y) : x 2 + y 2 < r 2 }, and by dD r its boundary {(x, y) : x 2 + y 2 = r 2 }. 
Let us consider the function V(x,y) = \ [x 2 +y 2 ). Its derivative along the 
solutions of ^ is 

V(x,y) = -y 2 <f)(x,y). 

Lemma 1. Let U be a bounded set, with a := sup{-\/a; 2 + y 2 , (x,y) G U}. If 
0(^2/) ^ out of U , and cj>(x,y) does not vanish identically on any dD r , for 
r > a, then every ^(t) definitely enters the disk D a and does not leave it. 

Proof. The level curves of V{x, y) are circumferences. For every r > a, 
the disk D r contains U. Since V(x,y) — —y 2 (j)(x,y) < on its boundary, 
such a disk is positively invariant. Let 7 be an orbit with a point j(t*) such 
that d* = dist (-f(t*),0) > a. Then 7 does not leave the disk D^*, hence it is 
positively bounded. Moreover j(t) cannot be definitely contained in dD r , for 
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any r > a, since V(x, y) does not vanish identically on any dD r , for r > a. 
Now, assume by absurd that j(t) does not intersect B a . Then its positive limit 
set is a cycle 7~(i), having no points in D a . The cycle 7~(t) cannot cross outwards 
any dD r , hence it has to be contained in dD r , for some r > a, contradicting 
the fact that V(x,y) does not vanish identically on any dD r , for r > a. Hence 
there exists t + > t* such that r y(t + ) G D a . Then, for every t > t + , one has 
-f(t) £ D a , because V(x, y) < on B a . $ 

Collecting the results of the above statements, we may state a theorem of 
existence and uniqueness for limit cycles of a class of second order equations. 
We say that an equilibrium point O is negatively asymptotically stable if it is 
asymptotically stable for the system obtained by reversing the time direction. 

Theorem 2. If the hypotheses of theorem^ and lemma^hold, and </>(0,0) < 0, 
then the system ^ has exactly one limit cycle, which attracts every non-constant 
solution. 



Proof. By the above lemma, all the solutions are definitely contained in D„. 
The condition </>(0, 0) < implies by continuity 4>{x,y) < in a neighbourhood 
No of O. This gives the negative asymptotic stability of O by Lasalle's invari- 
ance principle [TT], since V(x,y) > in N Q , and the set {V(x, y) = 0} n N Q — 
{y = 0} l~l No does not contain any positive semi-orbit. The system has just 
one critical point at O, hence by Poincare-Bendixson theorem there exist a limit 
cycle. By theorem [T] such a limit cycle is unique. £ 

This proves that every non-constant solution to the equation ^ studied in 
[S] is attracted to the unique limit cycle. 

We can produce more complex systems with such a property. Let us set 

n 

4>{x, y) = -M + ^2 H 2k (x,y), 

k=l 

with H 2 k(x, y) is a homogeneous function of degree 2k, positive except at O, M 
is a positive constant. Then, by Euler's identity, one has 

k=l ^ y ' k=0 

If <p(x,y) does not vanish identically on any dD r , for instance if H2k(x,y) — 
(x 2 + xy + y 2 ) k , then the corresponding system ^ has a unique limit cycle. In 
general, it is not necessary to assume the positiveness of all of the homogeneous 
functions H2k(x, y), as the following example shows. Let us set Q{x, y) = x 2 + 
xy + y 2 . Then take 

(f>(x,y) = + Q -Q 2 + Q 3 . 

One has 

v = x4> x + ycf> y = 2Q- 4Q 2 + 6Q 3 = Q(2 - AQ + 6Q 2 ). 
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The discriminant of the quadratic polynomial 2 — AQ + 6Q 2 is —32 < hence 
v > everywhere but at O. Moreover, <p(x, y) does not vanish identically on 
any circumference, hence the corresponding system Mm has a unique limit cycle. 



3 Non-linear g(x) 



Even if the equation we consider in this section are of a more general type, 
we actually derive our result from that of the previous section, so that we can 
consider what follows a corollary of the previous result. Let us consider the 
equation 

x + x$(x,x) + g(x) = 0. (12) 

We assume that xg(x) > for x ^ 0, g e C 1 (IR, H), g'(0) ^ 0, We could con- 
sider equations defined in smaller subset of the plane, without essential changes. 
The main tools is the so-called Conti-Filippov transformation, which acts on the 
equivalent system 

x = y, y = -g(x) - y$(x, y), (13) 

in such a way to take the conservative part of the vector field into a linear one. 
Let us set G(x) = J g(s)ds, and denote by cr(x) the sign function, whose value 
is —1 for x < 0, at 0, 1 for x > 0. Let us define the function a : 1R — ?► 1R as 
follows: 

a(x) = a(x)y/2G{x). 
Then Conti-Filippov transformation is the following one, 

(u,v) = A(x,y) = (a(x),y). (14) 

Since we assume that g £ C 1 (H, K) , one has a £ C 1 (H, H) . The function 
u = a(x) is invertible, due to the condition xg(x) > 0. Let us call x = /3(u) its 
inverse. The condition g'(0) > guarantees the differentiability of f3(u) at O. 
For x 0, that is for 11 ^ 0, one has, 

y/2G(x) a'{P{u)) g{x) 9{p{u)) 

For x = u = one has, 

a'(0) = Vff 7 (0) 5 /?'(0) ' 

Finally, 



<?'(0)' 



a^O u 
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Theorem 3. Assume g 6 C 1 (H, H), with g'(x) > and xg(x) > for x ^ 0. 
Let <t> e C^JR^IR 2 ) safe/?/ 



ct(x)v / 2G(^) 



2G(x) -r^ +$(x,y) 



TTien (751) /tas ai most one limit cycle. 



Proof. For u ^ 0, the transformed system has the form 



it = 



v = - 5 ( ( a( u ))- u $( i 9(t t ),u). 



(16) 



For it = 0, the above form is extended by continuity. In the following we consider 
only the case u ^ 0, that is x ^ 0, since the case u = = x is obtained by 



having the same orbits as ( 16 ) 



continuity. We may multiply the system ( 16 1 by g ^ u ^ > obtaining a new system 



v = —u — v- 



<?(/?(«)) 



(17) 



Such a system is of the type ([9|, so that we may apply theorem [TJ to get unique- 
ness of solutions. This reduces to require the strict star-shapedness of the func- 
tion 



that is, 



+ v$ v (P(u),v) > 0. 



Since /3(u) — x, v = y, the second term in the above sum is just y<S> y (x, y). As 
for the the first one, one has 



u^ u (f3(u),v)(3'(u) + t>(f3(u),v))g(f3(u)) - u*(/3(tt), v)g'(p(u))p>(u) 



Replacing j3(u) with x and applying the formulae (15) one has 



a{x)j2G{x) ( v CT (- T )V2GM^(a;,2/)Jg(.T) - a(x)^2G(x)^(x,y)g'(x) 



9{x) 
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Since a(x) 2 — 1 eveywhere but at 0, the above formula reduces to 



u$(J3(u),v) 

gWW) 

Concluding, one has 



2G(x) <S> x (x,y)g(x) - ${x,y)g'{x) ${x,y) 



9(x) 



20(a) vMx)-*{x, V)9'{x) + ^ y) 



9{xf 



Hence, the star-shapedness conditions reduces to 



a(x) v / 2G(x) 
9{x) 



2G(x f^y)9(x)-Hx,y)g'(x) + 
g{xy 



y$ y (x,y) ^ 0. 



x 2 a(x)j2G(x) 

If g(x) — x, then G(x) = — and = 1, for i ^ 0. In this case 

2 9\ x ) 

the star-shapedness condition just reduces to what considered in the previous 
section, since 

ont \^(x,y)g(x)-^(x,y)g'(x) 

2G{x) —-^ h $(ar,j/) = x$ x {x,y). 

g(xy 



Now we prove the non-linear analogous of lemma [T] Let us set 

E(x,y) = G(x) + ^- 

For r > 0, we set A r = {(x,y) : 2E(x,y) < r 2 } and dA r = {(x,y) : 2E(x,y) = 
r 2 } 

Lemma 2. Let U be a bounded set, with a := swp{y/2E(x, y), (x,y) € U}. If 
&(x,y) > out ofU, and &(x,y) does not vanish identically on any dA r , for 
r > a, then every j(t) definitely enters the set A a and does not leave it. 



Proof. Performing Conti-Filippov transformation, the sets A r are taken into 
the sets D r , as well as the boundaries dA r are taken into the boundaries dD r . 
The function <5>(x,y) does not vanish identically on any dA r if and only if the 
function <p(u, v) — $(/3(u), v) does not vanish identically on any dD r . Then one 
can apply lemma [I] to the system (17). In fact, the deri vative of the Liapunov 
function V(u, v) — \{u 2 + v 2 ) along the solutions of (17 1 is just 



V(u,v) 



i U$(/3(u),t>) 
9(P(u)) ' 



The function g ^ u ^ is positive for it 7^ 0, hence the hypotheses of lemmaJTjare 
satisfied by the system (17). As a consequence, the conclusions of lemmafmiold 
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for the system (17), and applyng the inverse transformation A 1 one obtains 
the thesis. A 

Now we can conclude proving the analogue of theorem [2] for the equation 
with a non-linear g(x). 

Theorem 4. If the hypotheses of theorem^ and lemma^ hold, and 0(0, 0) < 
0, then the system (13) has exactly one limit cycle, which attracts every non- 
constant solution. 

Proof. As the proof of theorem [2j replacing the Liapunov function V(x,y) 
with the Liapunov function E(x, y). 4k 
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